function analyze_results(theta2_range, rho_rho_star)
    % 分析散射结果
    
    % 找到主散射峰
    [peak_value, peak_idx] = max(rho_rho_star);
    peak_angle = rad2deg(theta2_range(peak_idx));
    
    % 计算半高宽（FWHM）
    half_max = peak_value/2;
    above_half = rho_rho_star > half_max;
    fwhm_indices = find(diff(above_half));
    if length(fwhm_indices) >= 2
        fwhm = rad2deg(theta2_range(fwhm_indices(end))) - ...
            rad2deg(theta2_range(fwhm_indices(1)));
    else
        fwhm = NaN;
    end
    
    % 显示分析结果
    fprintf('\n散射分析结果:\n');
    fprintf('主散射峰角度: %.2f°\n', peak_angle);
    fprintf('散射峰强度: %.2e\n', peak_value);
    fprintf('半高宽: %.2f°\n', fwhm);
end 